function[rfit,order]=gapopfitness(gaeachfitness,phono,limitted)
n=length(phono(:,1));
for i=1:n
    if limitted==0
        for j=1:length(phono(i,:))/2
            x(j)=phono(i,2*j)*10^phono(i,2*j-1);
        end
    else
        x=phono(i,:);
    end
    fit(i)=gaeachfitness(x);
end
[rfit,order]=sort(fit);
    